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ABSTRACT 

We report the detection of short period variations in the stars HD 69013 and HD 96237. 
These stars possess large overabundances of rare earth elements and global magnetic 
fields, thus belong to the class of chemically peculiar Ap stars of the main sequence. 
Pulsations were found from analysis of high time resolution spectra obtained with 
the ESO Very Large Telescope using a cross correlation method for wide spectral 
bands, from lines belonging to rare earth elements and from the Ha core. Pulsation 
amplitudes reach more than 200 ms^^ for some lines in HD 69013 with a period of 
11.4 min and about 100 ms~^ in HD 96237 with periods near 13.6 min. The pulsations 
have also been detected in photometric observations obtained at the South African 
Astronomical Observatory. 
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1 INTRODUCTION 

The rapidly oscillating Ap (roAp) stars are magnetic main 
sequence stars that pulsate in high radial overtone p modes 
with periods in the range of 5.7 — 21 min. They show 
broad-band photometric amplitudes less than 0.01 mag, 
whereas rapid radial velocity variations in rare earth ele- 
ment l ines can reac h several kms~^ (e.g.. lFrevhammer et al.l 
l2009l : iKurtj Il990l '). The roAp stars are important tar- 
gets for the study of the interactions among chemical 
anomalies, magnetic field and pulsations. These stars show 
abnor mal atmospheric st r ucture with chemical str atifica- 
tion (|Cowlev et al. I I2OOII : iRvabchikova et all I2OO2D . The 
pulsations of these stars also make them promising ob- 
j ects for interior model testing using ast eroseismology 
(|Aerts. Christensen-Dalsgaard fc Kurt j bOlOf ). The roAp 
stars were discovered by Kurt j (|l982l ): at present more than 
40 such stars are known. 

Interesting and surprising discoveries have been made in 
recent years, which give new insight in the study of pulsating 
Ap stars. New ground is being broken with /^ mag precision 
photometric data from the Kepler Mission. iBalona et al.l 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, as part of programmes 080.D-0191(A), 078. D- 
0192(A) and 078.D-0080(A) 



(|2010h discovered a roAp star that pulsates in both high 
overtone p modes and a low- frequency g mode, opening up 
the possibility of better modelling of the interiors of these 
most peculiar stars. 

In 2007 we started a high resolution survey of cool chem- 
ically p eculiar stars bas ed mostly on the photometric cata- 
logue o f lMartinej(| 19931 ). One of the goals of the survey is to 
select stars with high peculiarity and strong magnetic field, 
and to determine their fundamental parameters to select the 
most promising candidates to be roAp stars. Nearly 400 stars 
have been observed and many of them are good roAp can- 
didates for further high time resolution spectroscopic and 
photometric observations. For several of these stars we have 
obtained such observations and here we present the discov- 
ery of pulsations for two objects, HD 69013 and HD 96237. 

Both stars were in the lis t of stars with magnet i- 
cally split spectral lines found bv iFrevhammer et al.l(|2008D . 
HD 69013 is a typical cool Ap star; HD 96237 is more impres- 
sive with significant spectral variability. The physical param- 
eters Tefi and log g indicate that HD 69013 and HD 96237 are 
both main sequence stars, situated in the HR diagram where 
the instability strip crosses the main sequence. For many Ap 
stars in this region of main sequence rapid oscillations have 
been detected, therefore both stars were promising objects 
for pulsation testing. 
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Table 1. A journal of observations of HD 69013. The columns 
give the Julian Date (JD) of the start of exposure, the observation 
time and the number of spectra or photometric measurements. 



JD 


observation 


exposures 


Telescope 




time (min) 


number 




2454482.546 


57 


34 


ESO UT2 VLT 


2454502.549 


57 


34 


ESO UT2 VLT 


2455226.448 


121 


362 


SAAO 1-m 


2455229.377 


84 


252 


SAAO 1-m 



Table 2. A journal of observations of HD 96237. The columns are 
the same as for Table 1. 



JD 


observation 


exposures 


Telescope 




time (min) 


number 




2454540.515 


57 


34 


ESO UT2 VLT 


2455225.501 


118 


354 


SAAO 1-m 


2455228.438 


100 


301 


SAAO 1-m 


2455229.486 


90 


541 


SAAO 1-m 


2455235.483 


58 


184 


SAAO 1-m 


2455329.226 


224 


1256 


SAAO 0.5-m 



2 OBSERVATIONS AND DATA REDUCTION 

High time resolution spectroscopic observations were car- 
ried out at the European Southern Observatory (ESO) us- 
ing the Ultraviolet and Visual Echelle Spectrograph (UVES) 
installed at Unit Telescope 2 (UT2) of the Very Large Tele- 
scope (VLT). For HD 69013, data were obtained during two 
high time resolution observing runs on 2008 January 17 and 
February 6. For each run 34 spectra were obtained with 
exposure times of 80 s and readout and overhead times of 
~ 21s, corresponding to a time resolution of ~101s. For 
HD 96237 we obtained 34 spectra on 2008 March 15 with 
the same exposure and readout times. The wavelength re- 
gion observed is AA4970 — 7010 A, with a small gap in the 
region around 6000 A caused by the space between the two 
CCDs. The average spectral resolution is about R — 10^. 
The CCD frames were processed using the UVES pipeline 
to extract and merge the echelle orders to ID spectra that 
were normalised to the continuum. 

We also obtained photometric observations of HD 69013 
and HD 96237 in January, February and May 2010. These 
observations were obtained at the South African Astronom- 
ical Observatory (SAAO) in Johnson B filter with the 1-m 
telescope and SAAO CCD and STE4 detector and with the 
modular photometer at the 0.5-m telescope. The reduction 
of photometric observations was done with ESO-MIDAS soft- 
ware and with software developed at SAAO. The lists of the 
observations are shown in Tables 1 and 2. 



3 PULSATION SEARCH AND ANALYSIS 

For roAp stars lines of rare earth elements show higher 
pulsation amplitudes, while the lines of other chemical 
species, including hght elements and iron peak elements, 
show much s maller pulsation amplitude, or show none at 
all (see, e.g.. iMalanushenko. Savanov fc Rvabchikwalll998l : 



iKurtz. Elkin fc Mathvll2007^ . This strange behaviour is ex- 
plained by stratification where rare earth elements con- 
centrate in the upper layers of the stellar atmosphere 
where oscillation amplitudes reach a maximum, while most 
other chemical elements tend to concentrate in deeper lay- 
ers where the pulsation amplitude is lower. Lines of iron 
peak elements in roAp s tars show very low pulsation am- 
plitude or none at all (iKochukhov fc Rvabchikoval 1200 ll . 
lElkin. Kurtz fc Mathvij|2008l ). 

To search for rapid radial velocity variability we per- 
formed cross-correlation of sections of the spectrum using 
ESO-MIDAS software. We also measured the central posi- 
tions for profiles of individual spectral lines by the centre of 
gravity method. Frequency analyses of radial velocity and 
photometric time series were performed using ESO-MIDAS's 
Time Series Analysis an d a discrete Fourier transform pro- 
gramme bv lKurtzi (| 19851 ). 

3.1 HD 69013 

Cross correlation for the spectral band 5000 — 5800 A with an 
average spectrum taken as a template shows obvious rapid 
oscillations for HD 69013, as can be seen in Fig.[T] For the 
two independent observing data sets that we obtained with 
the ESO VLT telescope we find in the amplitude spectra 
highest peaks aX v = 1.46 mHz and v = 1.42 mHz, corre- 
spondingly, with a full-width-half-maximum uncertainty of 
0.07 mHz. Hence these two independent peaks are at the 
same frequency within the frequency error. 

The spectral lines of the rare earth elements also show 
pulsation with different amplitudes. The highest amplitude 
we detected was obtained for lines of Prill, shown in Fig.[5J 
while lines of Ndlll reveal a smaller amplitude as shown 
in Fig.O The other lines which belong to Eull, Cell, Lall 
also show pulsations with significant peaks in the amplitude 
spectra. 

The pulsation amplitude is smaller for the second ob- 
serving run, which suggests either rotational modulation 
or multiperiodicity. While many Ap stars show rotational 
light variations caused by abundance spots usually associ- 
ated with their magnet ic poles, there is no suc h evidence 
yet found for HD 69013 (|Frevhammer et al.|[200l ). This con- 
straint is weak and does not rule out rotational modulation 
pulsation amplitudes. Further observations are needed to 
study this question. 

The photometric observations of this star obtained at 
SAAO also show rapid oscillations, as can be seen in Fig.Hl 
Previo us photometric observations by iNelson fc Kreidj 
(|l993l ) and [Martinez fc Kurt j (|l994h did not detect varia- 
tions in this star. The amplitude spectra of our SAAO pho- 
tometric observations are shown in Fig. 4. We obtained data 
on two nights. The first data sets shows a signal at the same 
frequency as for the radial velocity data, but there is no 
significant peak in the second data set, as can be seen in 
Fig. 4. 

3.2 HD 96237 

The case for pulsation in HD 96237 is not as strong as for 
HD 69013, but all of the evidence together gives confidence 
that pulsation has been detected and this star is a roAp 
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Figure 1. Top: The amplitude spectrum of HD 69013 obtained 
from cross correlation over the spectral region 5000 - 5800 A for 
the first observing run (see table 1). An average spectrum was 
used as the template for the cross correlation. The frequency peak 
is at 1/ = 1.46 mHz. Bottom: the amplitude spectrum for the 
second spectroscopic run. The frequency peak \s u = 1.42 mHz. 
The pulsation amplitude dropped significantly between the two 
runs, indicating rotational modulation or multiperiodicity. 



star. The five panels of Fig. [5] show amplitude spectra for 
spectroscopic analysis of HD 96237. Clear pulsation peaks 
were obtained from cross correlation for the spectral region 
5000-5600 A using an average spectrum as a reference tem- 
plate, and for the core of the Ha line. 

As mentioned above, the spectrum of HD 96237 is very 
rich in rare earth element lines. We tried to measure pul- 
sations for individual spectral lines which were identified 
by comparison with synthetic s pectra calculated with the 
SYNTH code of lPiskunovl (|l992l '). but found that the pulsa- 
tion amplitude is too low for clear detection in most individ- 
ual lines, given the noise level in our spectra and relatively 
short observing run. In a few spectral lines of rare earth 
elements a peak corresponding the pulsation frequency de- 
tected by cross correlation is visible in the amplitude spec- 
trum, as can be seen in the third panel of Fig. [S] for a single 
Cell line. 

The combination of several good spectral lines reduces 
the noise level and produces a more reliable picture. The two 
lower panels in Fig. [S] show amplitude spectra for combina- 
tions of five rare earth elements lines giving a clear peak. 
For twelve lines of Fe I and Fe ll no pulsation is detectable, 
as is typical of the roAp stars. 

As a further test we examined a combination of several 



Figure 2. Top: the amplitude spectrum of HD 69013 obtained 
from the line of Pr III at A5299.993 A. The frequency peak is 
V = 1.51 mHz. Bottom: The amplitude spectrum for the second 
spectroscopic set of HD 69013 for the same line. The frequency 
peak is 1^ = 1.48 mHz. 



telluric lines and found that there is no signal above a noise 
level of 12ms~^. We conclude that the case that HD 96237 
is a roAp star is good. Further studies - particularly at the 
rotation phase where the pulsation amplitude is highest - 
will give more confidence. 

Photometric observations were obtained at SAAO for 
additional testing of pulsations in HD 96237. The frequency 
analysis of the photometric data listed in Table 1 is pre- 
sented in Fig.|Sl The upper panel of this figure supports the 
pulsation period found by spectroscopy , but other photo- 
metric observations sho wn in the bottom four panels and by 
iNelson fc Kreidj (| 19931 ) do not detect pulsation in this star. 

This can be understood in terms of rotational modu- 
lation in an oblique pulsator, and/or beating of multiple 
frequencies. The pulsation amplitudes in many roAp stars 
are modulated with the rotation period of the star. Thus 
some of the photometric observations may not have been 
in best aspect. To judge this the rotational period and 
ephemeris of HD 96237 needs to be determined. From the All 
Sky Automated Survey ( ASAS) and Hipparcos photometry, 
iFrevhammer et al] (|2008h found photometric variations with 
a period of 20.91 d, which may be the rotational period. Us- 
ing two seasons of dat a from the WASP (W ide Angle Search 
for Planets) Project (|Pollacco et al.ll2006h covering the in- 
tervals 2007 January 4 to 2007 June 3 and 2008 January 5 
to 2008 May 28 we obtained a similar period. Combining the 
WASP (passband from 400 to 700 nm) and ASAS (V-band) 
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Figure 3. Top: The amplitude spectrum of HD 69013 obtained 
from three lines of Ndm (AA5050.695 A, 5102.427 A, 5294.113 A) 
combined. The frequency peak is u = 1.47 mHz. Bottom; The 
amplitude spectrum for the second spectroscopic set of HD 69013 
for the same three lines. The frequency peak is = 1.44 mHz. 
Within the errors these two peaks are the same. 

photometry we find the following rotational ephemeris for 
the photometric maximum: 

HJD{phot.max.) = 2454483.9285 + 20.91S (1) 

The rotational phases calculated from this ephemeris 
are also shown in Fig.[Sl The bottom panel is for data ob- 
served in the same rotation period as the upper panel, but 
does not display any pulsations. The spectral observations 
presented at Fig. [S] were also obtained at a similar rotation 
phase using the above ephemeris. Assuming our case for pul- 
sation in this star to be good, we suggest two possible solu- 
tions: 1) the rotation period may be not 20.91 d, but double 
that. This ambiguity happen s for some peculiar stars (see 
for example lWade et al.l (|l997h '). Longitudinal magnetic field 
measurements over the rotational period can resolve this; 2) 
the star may be multiperiodic. There is some hint in Fig.[6] 
of a peak at 0.79 mHz. Additional observations are required 
to resolve these questions. 



4 CONCLUDING REMARKS 

The chemically peculiar magnetic star HD 96237 is an in- 
teresting object which has a very peculiar spectrum with 
significant spectral variability. The star demonstrates large 
overabundances of rare earth elements. A high resolu- 
tion spectrum obtained with the ESO 2.2-m telescope and 
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Figure 4. Top: The amplitude spectrum of HD 69013 obtained 
for the first run of SAAO photometry through a B filter. The 
highest peak has a frequency of 1.39 mHz. Bottom: The second 
photometric set for the star does not show any significant peak 
for frequencies close to that detected by spectroscopy. 

FEROS spectrograph resembles the spectrum of another 
highly peculiar star, HD 101065. Abundances of Ndll and 
Ndm determined from this spectrum are even higher than 
in HD 101065. Other rare earth elements also show large 
overabundances similar to those found in HD 101065. An- 
other spectrum of the star obtained with VLT UVES 
was significantly different from the FEROS spectrum with 
much less intense spectr al Unes of the rare earth elements 
(|Frevhammer et al.ll2008t ). The similarity of the spectra and 
chemical abundances of HD 96237 to those of HD 101065, 
which was the first detected rapidly oscillating Ap star, in- 
creases the interest of this star. It is not known what de- 
termines the pulsation amplitude of the roAp stars. While 
the principal periods of HD 96237 and HD 101065 are similar 
(13.6 min and 12.1 min, respectively), their photometric and 
spectroscopic radial velocity amplitudes are very different. 

Our results for HD 69013 demonstrate clearly that this 
is a new roAp star with a low pulsation amplitude. This 
star shows possible rotation modulation and may be a useful 
target to study pulsation behaviour over rotation period. 
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Figure 5. A cross correlation and Fourier analysis for 34 spectra 
of HD 96237. From top to bottom: the first panel shows a highest 
peak with a frequency of 1.22 mHz; the second panel has a peak 
at 1.16 mHz which is in consistent with the first panel within the 
errors. Both upper panels also show peaks with the frequency 
0.78 mHz which is close to that detected in two photometric ob- 
servations (see panels 3 and 4 at Fig.|6] The peaks in the third 
and forth panels have highest peaks at a frequency of 1.27 mHz. 
The fourth panel shows an amplitude spectrum for a combination 
of lines - Nd III 5102.43 A, 5294.11 A, Prill 5284.69 A, 5299.99 A, 
Ndll 5319.81 A. No pulsation signal was detected for a combina- 
tion of 12 lines of iron, as shown at the fifth panel. 
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Figure 6. Amplitude spectra of the photometric observations 
for of HD 96237. Extinction and sky transparency variations were 
removed using a second degree polynomial. The data are ordered 
in time from top to bottom. The highest peak in the upper panel 
has amplitude near 4mmag and frequency of 1.19rhHz, which is 
at the error limit of the frequency found for the spectroscopic 
data. The star may also show a harmonic with a frequency of 
2.42 mHz. The third and fourth panels show peaks with the same 
frequency of 0.79 mHz. A similar frequency peak was detected in 
spectroscopic observations (see Fig.[5]l. It is not clear weather this 
peak is real; only one other roAp star, HD 116114, with a s imilar 
low frequency has previously been found l lElkin et al.ll2005lV 



Strasbourg, France. This paper uses observations made at 
the South African Astronomical Observatory (SAAO). 
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